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ABSTRACT 

Nuclear and Coulomb collisions in GRB jets create a hot e plasma. This collisional heating 
starts when the jet is still opaque and extends to the transparent region. The e ± plasma radiates 
its energy. As a result, a large fraction of the jet energy is converted to escaping radiation with 
a well-defined spectrum. The process is simulated in detail using the known rates of collisions 
and accurate calculations of radiative transfer in the expanding jet. The result reproduces 
the spectra of observed GRBs that typically peak near 1 MeV and extend to much higher 
energies with a photon index j3 ~ —2.5. This suggests that collisional heating may be the 
main mechanism for GRB emission. 

Key words: gamma-rays: bursts, theory — plasmas — radiation mechanisms: thermal, non- 
thermal — radiative transfer — relativity — scattering. 
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1 INTRODUCTION 

Cosmological gamma- ray bursts (GRBs) are associated with ultra- 
relativistic jets from short-lived powerful sources such as hyper- 
accreting, just-born black holes. The jet starts as an opaque black- 
body fireball that accelerates, expands and releases its thermal ra- 
diation at the photospheric radius R*. One may expect a quasi- 
blackbody spectrum from such jets (Paczynski 1986; Goodman 
1986), similar to the relict radiation from the big bang. However, 
the simpe blackbody model is inconsistent with observations (e.g. 
Preece et al. 2000). It is clear that some form of heating operates 
in the jet and changes its radiation from blackbody. Heating may 
occur at radii r < and change the photospheric radiation via 
Comptonization. It may also occur at radii r > J?* and generate 
nonthermal synchrotron emission. 

Two heating mechanisms are usually considered in GRB jets: 
internal shocks and dissipation of magnetic energy. The details of 
both mechanisms are uncertain as they depend on complicated col- 
Usionless processes in the plasma. A long-standing problem is the 
radiative efficiency of these processes. 

Recent observations by Fermi telescope provided new data in 
a broad spectral range from 8 keV to ~ 100 GeV. The data confirm 
the previous BATSE result that the prompt GRB spectrum typi- 
cally peaks near MeV (Preece et al. 2000). The typical spectrum 
is approximately described by the Band function (a smoothly bro- 
ken power law), which extends to high-energy bands with a photon 
index f) ~ —2.5. The prompt GRB radiation is highly variable 
on timescales as short as millisecond, suggesting a small radius of 
emission, possibly comparable to the photospheric radius 



The spectrum of GRB emission at very high energies E J; 
1 GeV remains so far uncertain. Multi-GeV photons overlapping 
the prompt MeV radiation have been detected in about 5 per cent 
of GRBs, however they may be produced by a distinct source at 
large radii, e.g. by the blast wave from the explosion. The distinct 
source is visible in several GeV-emitting bursts (e.g. Abdo et al. 
2009b; Ryde et al. 2010) and likely present in all of them, obscuring 
the behavior of the prompt Band spectrum at high energies (cf. the 
debate over GRB 080916C: Abdo et al. 2009a; Kumar & Barniol 
Duran 2009; Ghisellini et al. 2009). 

The standard theory (e.g. Paczynski 1990) predicts that jets 
with r ~ 10 3 must produce bright photospheric emission. Ob- 
servational search for this emission usually assumed that it has a 
quasi-thermal spectrum (e.g. Ryde 2005). It was argued that the 
absence of the photospheric emission component would imply that 
the jet is magnetically dominated and cold, with negligible initial 
thermal energy (e.g. Daigne & Mochkovitch 2002; Zhang & Pe'er 
2009). 

The main finding of the present paper is that the Band-type 
spectrum naturally forms in GRB jets as a result of collisional heat- 
ing. Most of radiation produced by this mechanism is emitted near 
the photosphere i?» and therefore it may be called photospheric to a 
first approximation. Our result supports the view that photospheric 
emission is not a rare quasi-thermal component; instead, it is the 
main Band component of GRB emission that is routinely observed 
in all bursts. 

This paper considers the standard model of a baryonic jet with 
comparable numbers of neutrons and protons (our model would 
not work for jets that are completely dominated by magnetic field, 



* Also at Astro-Space Center of Lebedev Physical Institute, Profsojuznaja 
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1 Recent suggestions that the prompt 7-ray emission must come from a 



large radius r S> R* (e.g. Racusin et al. 2008; Abdo et al. 2009a) are based 
on incorrect assumptions (see Section 6.1). 
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with negligible baryonic loading). Before becoming transparent to 
radiation, the jet evolves to the two-fluid or 'compound' state: a 
plasma with bulk Lorentz factor T embeds a neutron flow with 
Lorentz factor r„ < T (Derishev, Kocharovsky & Kocharovsky 
1999a; Bahcall & Meszaros 2000; Fuller, Pruet & Abazajian 2000; 
Meszaros & Rees 2000; Rossi, Beloborodov & Rees 2006; Koers 
& Giannios 2007). Regardless the details of their formation, com- 
pound jets with r„ < F have a robust feature: nuclear collisions 
between the neutron and proton fluids continually create multiple 
e with energies ~ m-n-c 2 ~ 140 MeV. Their energy is immedi- 
ately converted to radiation before e join the thermalized plasma. 
Nuclear collisions also heat the proton component of the jet, and 
protons gradually drain their energy into thermalized e plasma 
via Coulomb collisions. 

Similarly to internal shocks, collisional heating taps the ki- 
netic energy of internal motions in the jet — the streaming of 
plasma through the neutron component with a relative Lorentz fac- 
tor r re i = |(r/r n + r n /r). In contrast to internal shocks, the 
heating is not confined to a shock front. It operates in volume. 

Several works previously proposed that some sort of volume 
heating shapes the spectrum of GRB emission (e.g. Thompson 
1994; Ghisellini & Celotti 1999; Stern & Poutanen 2004; Rees & 
Meszaros 2005; Pe'er, Meszaros & Rees 2005; Giannios & Spruit 
2007; Ioka et al. 2007; Asano & Terasawa 2009). The models as- 
sumed some form of collisionless dissipation, which is difficult to 
calculate from first principles. In this context, two special features 
of our model should be noted: 

(i) The collisional heating is robust, and its history in the expand- 
ing jet is well defined. The rate of collisions determines the radial 
dependence of the heating rate Q oc r~ 2 . 

(ii) The collisional heating injects energy into e via two branches 
with comparable heating rates: 

(a) Nuclear collisions maintain a continual e cascade in the 

(b) Coulomb collisions in the two-temperature plasm^l con- 
tinually transfer energy from protons to thermalized e . 

Branch (a) is important because it loads the jet with a large number 
of e pairs and determines the photospheric radius of the burst. On 
the other hand, it will be shown that branch (b) plays an important 
role in the formation of the GRB spectrum. The radiation emerging 
from a collisionally heated jet has a well defined spectrum, which 
can be calculated numerically. This radiative transfer problem is 
solved in this paper using a Monte-Carlo code that tracks the evo- 
lution of photons and e in the heated and expanding plasma flow. 

The paper is organized as follows. Section 2 gives a compact 
summary of the simplest model of GRB jets with no internal dissi- 
pation. Such jets passively cool down as they expand, and eventu- 
ally emit thermal radiation whose spectrum cuts off exponentially 
at ~ MeV. It is inconsistent with the observed GRB spectra. We use 
the model of a passively cooling jet as a benchmark and a first test 
problem for our radiative transfer code. 

Section 3 describes neutron-loaded jets and formation of com- 
pound flows with r n < r. Section 4 describes the collisional radia- 
tive mechanism operating in compound flows. Section 5 presents 
the radiation spectrum received by distant observers. The results 
are discussed in Section 6. Section 6 also discusses the possibility 
of additional emission that may be generated by neutron decay. 



jet. 



GRB outflows are believed to be beamed and therefore called 
'jets' throughout this paper. However, the results apply equally well 
to spherically symmetric outflows. As long as the opening angle of 
the explosion exceeds 1/T, the jet near the axis is causally dis- 
connected from its edge, and its dynamics is the same as that of a 
spherically symmetric flow. 



2 THERMAL EMISSION FROM PASSIVELY COOLING 
JETS 

We focus in this paper on jets that are accelerated by thermal (radi- 
ation) pressure, with a subdominant magnetic field. This standard 
model is briefly summarized below (see e.g. Paczyhski 1990). 

At small radii r, the GRB jet is in thermodynamic equi- 
librium and its luminosity is carried mainly by radiation L ~ 
(4/3)caT 4 r 2 4?rr 2 (hereafter we use the isotropic equivalent of lu- 
minosity, which would be produced by a spherically symmetric out- 
flow of the same density and temperature). As the jet expands adi- 
abatically, the ratio of photon and baryon number densities n 7 /n 
remains constant, i.e. effectively the photon number is conserved 
(similar to the cosmological big bang). The jet accelerates until the 
radiation energy density U~, = aT 4 decreases below the rest-mass 
density nm p c 2 . Then its Lorentz factor saturates at the asymptotic 
r. We will denote the characteristic saturation radius by R s . For a 
radially expanding jet, 



R s ~ Tro, 



(1) 



where ro is the radius at the base of the jet, at the beginning of its 
acceleration. The photon-to-baryon ratio in the jet is given by 



240rr^L- 

lb 

and the jet energy per photon is 



- TmpC n 
Eo ~ 



4,-^MeV. 



(2) 



(3) 



If no heat is generated by any dissipative processes, the ther- 
mal radiation trapped in the opaque flow continues to cool adia- 
batically at r > R 3 until it is released at the photosphere i?*. 
Between R B and i?*, the radiation temperature decreases as n* 1-1 
where -y = 4/3 is the adiabatic index of radiation, which gives 
E oc r~ 2 / 3 . The plasma has a small heat capacity (for n <C n 7 ) 
and simply tracks the radiation temperature. Electrons are ther- 
mally coupled to radiation via Compton scattering, and ions main- 
tain thermal equilibrium with electrons via Coulomb collisions at a 
common low temperature kT S3 1 keV in the rest frame of the jet. 

The optical depth of the jet is given by 



tt = 



4-7rrm p c 3 r 3 



n 1 £52 r 



52 1 3 , 



(4) 



2 The thermalized are Compton-cooled and kept at a temperature much 
lower than the proton temperature (see Section 4.4). 



where or = 6.65 x 10~ 25 cm 2 is Thomson cross section and L = 
4nr 2 r 2 nm p c 3 is the isotropic equivalent of the kinetic luminosity 
of the jet (it approximately equals the total luminosity at r > R s )- 
The photosphere radius R* ~ 10 10 L52r^" 3 cm is larger than R s 
for T < 10 3 L^rg^' ' 4 . The radiation luminosity released at the 
photosphere is L 7 ~ (R+ / R s )~ 2 ^ L, and the mean energy of the 
escaping photons is given by 

E(IL) ~ " 2/3 E « 4Fr 3 L 5 - 2 5/1 V^ 7 6 MeV. (5) 

The thermal radiation creates a bright burst with MeV peak if T ~ 
10 3 . The burst is weak for slower flows: L 7 oc T 8 ^ 3 . 
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Figure 1. Photon spectrum emitted by a passively cooling jet to distant 
observers. The jet has kinetic luminosity L = 10 52 erg/s (isotropic equiv- 
alent) and initial size ro = 10 7 cm. It implies the total number flux of 
photons (isotropic equivalent) TV = 1.6 X 10 57 s — 1 . The photon spec- 
trum has been calculated using the Monte-Carlo radiative transfer code. 
Two cases are shown: jets with asymptotic Lorentz factors V = 600 and 
300. The slopes of the spectra near 10 keV are close to 1.4, which corre- 
sponds to photon index a = 0.4 The spectrum cuts off exponentially at 
E > B pcak ~ 1 MeV. 



We have calculated the radiation spectrum emerging from the 
passively cooling jet using the Monte-Carlo code described in Ap- 
pendix B. The scattering of the initial Planck spectrum was fol- 
lowed until the jet expanded to complete transparency. The spec- 
trum of radiation received by distant observers is shown in Figure 1 
for two example models. In qualitative agreement with the above 
estimates, the jet with T = 600 produces a bright GRB, whose 
spectrum peaks near MeV and cuts off exponentially. Two details 
are worth mentioning: 

(i) Radiation emitted by passively cooling jets is not Planck- 
ian. Observer sees different parts of the photosphere at different 
angles, with different Doppler shifts. As a result, the low-energy 
slope of observed spectrum is softer than the Rayleigh- Jeans a = 1 
(photon index). Isotropic emission in the jet frame would produce 
a — 0, and the exact a is controlled by the photon angular distri- 
bution near photosphere (e.g. a w 0.4 near 10 keV in Fig. 1). 

(ii) The standard description of adiabatic cooling predicts that 
radiation from a source at optical depth tt is cooled by the factor 
A(tt) = r^ 2 ^ 3 by the time the jet expands to its photospheric ra- 
dius. The detailed transfer calculations give larger A. For example, 
A(8) = 0.58 instead of 8 _2/3 = 0.25 and A(20) = 0.39 instead 
of 20" 2/3 « 0.14. The scaling A oc t t 2,/3 is maintained at large 
tt > 10. 

The radiative transfer simulations illustrate and refine the stan- 
dard fireball picture. They show that thermal radiation emitted by a 
passively cooling jet with T > 500 peaks at 73 pC ak ~ 1 MeV and 
carries away a significant fraction e of the jet energy (e.g. e w 1/4 
for the model with T = 600 in Fig. 1). However, its spectrum 
cuts off exponentially at E > 73 pC ak. Therefore, the model fails 



to explain the observed GRBs, whose spectra extend to energies 

E -Epcak- 



3 NEUTRON COMPONENT AND COMPOUND FLOWS 

The picture described in Section 2 is incomplete, because it ne- 
glects the neutron component of the jet. In any plausible scenario 
of the GRB trigger, the central engine is dense, hot and neutron 
rich (Derishev, Kocharovsky & Kocharovsky 1999b; Beloborodov 
2003). In particular, accretion-disc models for GRBs predict a high 
neutron fraction (see Beloborodov 2008 for a review). The neutron- 
rich matter is expected to enter the relativistic jet (although the de- 
tails of this process may vary, see e.g. Levinson & Eichler 2003; 
Metzger, Thompson & Quataert 2008), and the neutron-proton jet 
initially accelerates as a single fluid where n and p are coupled by 
frequent nuclear collisions. 

Neutrons and protons tend to combine into helium where the 
jet temperature drops to 140 keV. This process, however, competes 
with rapid expansion and is marginally successful. Collimation of 
the jet generally helps nucleosynthesis because it slows down ex- 
pansion (cf. Fig. 4 and 5 in Beloborodov 2003). However, even in 
cases where helium production is complete, some neutrons survive 
in jets with a neutron excess, as helium production consumes equal 
numbers of n and p. In addition, free neutrons are produced by 
spallation of a-particles at larger radii where the jet is heated. 

The expansion of neutron-loaded jets generally leads to the 
formation of a compound flow, a slower neutron component with 
Lorentz factor T n is embedded in a faster proton flow with Lorentz 
factor r > T n . The compound flow develops at the characteristic 
radius R n where the timescale for n-p collisions becomes longer 
than the jet expansion time. 

In particular, in jets that accelerate to T > r cr it « 
400Z/52 4 r ( ^'' 4 neutrons do not develop the full Lorentz factor — 
their Y n saturates at a smaller value (Derishev et al. 1999a; Fuller et 
al. 2000). For example, a baryonic flow accelerated to T = 10 3 can 
contain neutrons with T n ~ 10 2 . In spite of the significant differ- 
ence in Lorentz factors, the two components move together without 
any separation for a long time, because their velocities relative to 
the central engine are almost equal (the velocities practically equal 
c). 

Even in jets with T < r cr it, compound flows with r n < V are 
expected to form, because the jet is variable. The neutron compo- 
nent does not participate in internal shocks that develop in variable 
jets. As a result, neutrons from the slow portions of the jet migrate 
across the shocks and penetrate the faster portions (Meszaros & 
Rees 2000). This internal mixing is caused by the short-timescale 
variability of the central engine that creates a non-uniform flow. 
The mixing occurs on scales 5r ~ r/T 2 , much smaller than the 
total thickness of the ejected flow. Large variations of Lorentz fac- 
produce a non-uniform compound flow with T/T n ^> 1, 
Neutron migration can be illustrated by the following simple 
model. Suppose neutron-loaded flow A with Lorentz factor Fa is 
followed by faster flow B with Lorentz factor Yb 2> Ta (Fig. 2). A 
shocked region C forms between the two flows; the shocked plasma 
has Lorentz factor Tc such that Ta < Tc < Tg. Initially, the neu- 
tron component of flow A is coupled to protons by frequent colli- 
sions, so they behave as a single fluid. When neutrons in flow A 

3 Large variations on small scales are suggested by observed variability in 
GRBs. Large variations also generally help explain the high efficiency of 
dissipation of internal motions in the jet (e.g. Beloborodov 2000). 
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rjTl : r/F c 2 

Figure 2. Faster flow B sweeps slower flow A and compresses it into a 
shocked shell C. Neutrons from flow A are not swept and instead penetrate 
flow B. As a result a compound flow is formed: flow B contains a slower 
neutron component with T n = Y A . The penetration depth of neutrons is 
~ r/rj, in the lab frame; it is (r^/r^) 2 larger than the thickness of 
shocked region C. 

become collisionally decoupled, they are not swept into region C 
anymore. Instead, they penetrate region B with the relative Lorentz 
factor r rc i = \{Tb/T a + Ta/Tb) » r B /2r A - The penetra- 
tion/mixing length is ~ (Tc /r^t) 2 larger than the thickness of the 
shocked region C. 

Some of the penetrating neutrons collide with their new host 
flow. Each collision dissipates the relative kinetic energy (r re i — 
l)m p c 2 . The number of collisions per baryon of flow B during 
the jet expansion timescale equals the collisional 'optical depth' of 
the slow neutrons r„ = n n ar/T„, where a is the nuclear cross 
section. At the beginning of neutron penetration r n ~ 1 and a large 
heat is generated by collisions. The collisions decelerate flow B 
from Tg to a new T, which is found from energy conservation in 
the static lab frame, r n r 2 /2r^ ~ TgQThis gives T that is lower 
than the original Tg by the factor (r n rs/2r^) -1 ' 2 as long as 

Tn>T-]. 

In summary, GRB jets are expected to contain a significant 
neutron component (unless they are essentially baryon-free and 
completely dominated by Poynting flux). At the radius R n where 
T n ~ 1, collisions between neutrons and protons become rare and 
compound flows with T > T n inevitably develop. The schematic 
picture in Figure 3 indicates the main characteristic radii of the jet. 
The rare nuclear collisions in the region r„ < 1 dissipate huge 
energy, comparable to the total energy of the jet. The dissipation 
efficiency of collisions is (r re i — 1)t„. It may exceed 100 per cent 
as the collisionally decelerated jet tends to regain its initial Lorentz 
factor via adiabatic cooling and re-dissipate its energy. Below we 
explore how collisional dissipation affects the jet radiation. 



4 RADIATIVE MECHANISM 

Hereafter we consider a simplified jet model: a neutron component 
with a single bulk Lorentz factor T n is embedded in a fast proton 




Figure 3. Schematic picture of a baryonic jet. The jet starts to accelerate at 
radius r'o . Compound flow with r n < T forms at radius R n (eq. |6) and 
strong collisional heating begins at this radius. The jet becomes transparent 
to radiation at the photosphere i?* ~ 20R n (eq. [26}; its position is regu- 
lated by creation in the heated region. The figure also shows the mean 
radius of neutron decay, Ra = 3 X 10 15 (r n /100) cm, and radius R^ ac 
where the jet starts to decelerate because of the interaction with an external 
medium. The photospheric emission is released at R t , and its spectrum is 
strongly modified by sub-photospheric collisional heating. Collisional heat- 
ing continues at r > R* , although with a smaller rate. 



component with constant Lorentz factor F 3> r„. The proper den- 
sities of the neutron and proton components will be denoted by n„ 
and n, respectively. 



4.1 Inelastic nuclear collisions 



We consider collisions at radii where r n = n n ar/r n < 1, 




where L n = / Fnr' 1 Y%ri n m v c z is the isotropic equivalent of the ki- 
netic luminosity of the neutron flow, and a~3x 10~ 26 cm 2 is the 
effective cross section for nuclear collisions. The rate of collisions 
per unit volume (a Lorentz-invariant quantity) is given by 



n — nn n r ro i(jc. (7) 

Here r re i = §(r/r„ + r n /T) « r/2r„ is the relative Lorentz 
factor of the neutron and proton components of the jet. 

Collisions between neutrons and protons occur with signifi- 
cant r rc i and hence have a large inelastic fraction f- me \ Jj 1/2 
(Amsler et al. 2008). The energy /i ne ir re im p c 2 converts to mildly 
relativistic pions. The data on ^ multiplicity in p-p collisions are 
found in e.g. Breakstone et al. (1984) and refs. therein; a similar 
multiplicity is expected for n-p collisions. The total 7T and ty° 
multiplicity is larger by the factor of 3/2; it is typically 5-6 for GRB 
jets. 

The pions immediately decay: ^ — > jx + — > e + Ve 
and 7T° — > 7 + 7- The produced neutrinos escape with observed en- 
ergies ~ O.ir GeV and carry away a fraction f u ~ 1/2 of the pion 
energy0This multi-GeV neutrino emission is an important predic- 
tion of the baryonic jet model (Derishev et al. 1999a; Bahcall & 
Meszaros 2000; Meszaros & Rees 2000), which may be verified 



4 The decelerated flow with T < Tg stores the heat of ~ (r/2r^)m p c 2 
per baryon, and later tends to regain its initial Lorentz factor Tg as the heat 
converts back to bulk kinetic energy via adiabatic cooling on the expansion 
timescale. 



5 On average, neutrinos take ~ 3/4 of n energy. The average fraction 
of 71-^ and 7r° energy that is given to neutrinos may be estimated as /„ ~ 
(2/3)(3/4) = 1/2. 
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or dismissed by future neutrino detectors. The existing upper lim- 
its from Super-Kamiokande experiment are ~ 10 times above the 
expected neutrino flux (Fukuda et al. 2002). 

The fraction 1 — /„ of pion energy is given to relativistic e 
and high-energy 7-rays, which quickly convert to e via 7-7 reac- 
tion. Thus, the net result of one inelastic collision is the injection of 
several e with a Lorentz factor 70 ~ m,r/m e « 300 in the rest 
frame of the plasma flow. The injected e carry a significant frac- 
tion f± — /inci(l — ft/) ~ 1/4 of the collision energy r rc im p c 2 . 

Since neutrons move with a negative radial momentum in the 
frame of the proton flow, e are injected with a negative momen- 
tum. However, they become quasi-isotropic in the plasma frame 
after one Larmor rotation in the magnetic field of the jet (the field 
is transverse to the jet direction; its radial component is strongly 
suppressed as follows from magnetic flux conservation in the ex- 
panding flow). Any reasonable magnetic field advected from the 
central engine implies a very short gyration timescale yom e c/eB, 
and the net momentum of injected e is immediately communi- 
cated to the plasma and vanishes in the plasma frame. 

The energy of e injected per unit volume per unit time equals 
/±r re im p c 2 h. Practically all of this energy is quickly converted 
to radiation (see Section 4.2). As a result, radiation accumulates 
energy density with rate (measured in the plasma rest frame), 



Qnth = /±r rc im p c h. 



(8) 



Here subscript 'nth' stands for radiation produced by nonthermal 
that are generated by nuclear collisions. Using equation {7} and 

Qnth = cT dQ nth /dr, one finds 

1 dQ n th f± r 



dlnr 



4 r„ 



16 FT 



Rn 



(9) 



is the dimensionless 'compactness' parameter of the radiation field. 
One can express I as 



I = e T„ 



U-, 



(13) 



Here e is the fraction of the jet energy that is carried by radiation 
(e > 0.1 for the model proposed in this paper), and 



OTiir 



H. JL En ~ ( —\ ( r 



GrO 



(14) 



The compactness / is high in the main heating region r n ^ 1, and 
hence Compton cooling is fast, tc <C t cxp . 

The high compactness has another implication. Photons that 
are scattered by relativistic e to energies E' S> 1 MeV in the jet 
frame will not survive - they will convert to secondary e via reac- 
tion7 + 7 — > e + + e~. The development of e cascade that accom- 
panies Compton cooling of relativistic particles (Appendix A) is 
described in detail by Svensson (1987) and Lightman & Zdziarski 
(1987). In this paper, the cascade is modeled numerically with our 
Monte-Carlo code. The typical multiplicity of secondary e M s 
is comparable to 60. The total multiplicity of created following 
one nuclear collision is 



M = MoMs ~ W , 



where 



Mo = 



/±r re im p 



3 T_ 

4 r n 



(15) 



(16) 



70 m e 

is the multiplicity of 'primary' e injected with Lorentz factor 
70 ~ m-n/me following a nuclear collision. 



4.2 Radiative cooling of injected e 

The pairs injected by nuclear collisions have a large Lorentz fac- 
tor 70 ~ m^/m e and immediately radiate their energy via Comp- 
ton and/or synchrotron cooling. 



4.2.1 Compton cooling 

The timescale for Compton cooling of an electron with Lorentz 
factor 7 by radiation with energy density U-y i jf| 



tc 



3m e c 



(10) 



4£/ 7 <tt7 

Radiation is initially present in GRB jets (Section 2), and (7 7 is 
further increased as the radiation absorbs the energy of injected 
e ± . Compton cooling timescale is shorter than the timescale of jet 
expansion t C x P = r/cT. Their ratio is 

tc 3 



tcxp 

where 



4Z 7 ' 



?7 7 
m e c 2 



(11) 



(12) 



6 This estimate assumes Thomson scattering, i.e. neglects the Klein- 
Nishina correction to the scattering cross section. The peak of GRB radi- 



ation is E' 



peak 



MeV/r ~ keV in the jet frame. Since fi' ^'y < m e c 



for all 7 70, most of the scattering by occurs in Thomson regime. 
Exact calculations of radiative transfer with the full Klein-Nishina cross 
section are performed in Section 5. 



4.2.2 Synchrotron cooling 

In the presence of magnetic fields, the injected e also experience 
synchrotron losses. The synchrotron cooling timescale is similar 
to equation \10\ except that [7 7 in this equation is replaced by the 
magnetic energy density Ub = B 2 /8n (measured in the jet frame). 
The synchrotron losses dominate if Ub > U-y. Ub may be ex- 
pressed as 



8tt 47rr 2 rV 



(17) 



where e B is the magnetic fraction of the jet energy. The typi- 
cal energy of synchrotron photons in the plasma frame is E' s — 
0.3 ^y 2 heB jm e c. The corresponding energy in the lab frame, E s ~ 
TE' S , is given by 



E s 



0.3 7 o 



h e 
m e cr 



1/2 



200r 12 1 4 /2l '52 2 k eV, (18) 



where we substituted 70 ~ m^/m e « 300. Synchrotron emission 
peaks in the region r J; R„ where heating peaks and most of e 
are injected. Jets for which synchrotron cooling is significant (i.e. 
where it can compete with Compton cooling) have large es\ then 
E s is comparable to the typical -Z5 pca k of observed GRB spectra. A 
similar E s ~ Speak was found by Koers & Giannios (2007). This 
feature of collisionally heated jets offers an additional mechanism 
for the preferential peak position at 0.1-1 MeV. 

Synchrotron emission from particles with low 7 5 is self- 
absorbed. These particles cannot be cooled by the synchrotron 
mechanism; they are Compton cooled. 



6 Andrei M. Beloborodov 



4.3 Optical depth of the jet 

In view of the strong e loading and the large cross section for 
photon scattering err ^> o, one may expect a large optical depth 
tt where the bulk of nuclear collisions occur. Note, however, that 
t t oc T" 3 while t„ oc T~ 3 . In compound flows (F/F„)~ 3 < 1; 
as a result, the collisionally heated jet with r n 1 has a moderate 
tt. 

If no pairs were created, the optical depth of the compound 
flow would equal t v (eq. IT4b . which may be smaller than unity at 
r ^ R n ■ The actual optical depth is enhanced and dominated by e ± 
created by the nonthermal cascade (Derishev et al. 1999a). The con- 
tinually injected e quickly cool down and accumulate at relatively 
low energies, forming a thermalized population that maintains a 
Maxwellian distribution via frequent Coulomb collisions between 
e . It is convenient to express the rate of e ± supply as 



n± = Mh ■■ 



YQ n 



(19) 



where Qnth = f±T Te im p c h is the rate of energy injection in pri- 
mary e , and Y = M s /"fo is the 'pair yield' of the cascade. A 
minimum Y m i n ~ would be obtained when counting only 
the primary e from pion decay. This may be appropriate for very 
strongly magnetized jets where synchrotron cooling of e strongly 
dominates over Compton cooling (Ub Uj) and suppresses the 
e ± cascade. In weakly magnetized jets with Ub < U-,, the devel- 
opment of e ± cascade gives Y = Aia/'Jo ~ 0.2. 

Let n± be the density of accumulated thermalized pairs. Their 
annihilation rate is given by 

Wann = ^ <T T Cn±. (20) 

16 

This expression assumes n± > n and kT < m e c 2 ; both assump- 
tions are valid where annihilation is significant. The density of ac- 
cumulated e evolves according to equation 

^-f- (r 2 n±) =n± — n an n- (21) 

At the beginning of collisional heating, h± and n ann are both larger 
than the left side of equation d21l l. and the equilibrium h± ~ n an n 
is established, 



Y J±i. rd acn„n = — <t t cn ± , 

m e 16 

which gives 

n±a T r /4 m p a T , , , L \ 1/2 F n 
T V 3 m e a L„ J Y 



(22) 



(23) 



/y\ 1 ' 2 / l \^ 2 1 r x- 1 

(02) {nj fe) - (24) 



or, using r n = R n /r and f± ?a 1/4 

-rr(r) = to — , to « 20 
r 

tt stays near the equilibrium value oc r n even after the annihilation 
timescale becomes long and the e population freezes-out. This 
is the result of a coincidence: the annihilation equilibrium gives 
tt oc r , which is also maintained when h± — n arm = 0. There- 
fore, equation \24\ remains valid even at late stages when the jet 
becomes transparent to radiation. 

The e optical depth is maximum at the beginning (and peak) 
of collisional dissipation when T n ~ 1. At this stage, r/r n is lim- 
ited by the deceleration effect of collisions (see the end of Sec- 
tion 3). In particular, the collision of flows A and B considered in 
Section 3 gives a compound flow with L ~ Lb, L n ~ La, and 
r/r„ ~ (r e /r^) 1/2 at r ^ R n . Then equation ([24} gives 



ro ~ 20 



Y_\ 1/2 (Ta Lb\ 1/2 



V Tb La J 



(25) 



It is reasonable to suppose Lb /La > 1 when Tb/Ya > 1 and 
expect (TbLa/YbLb) 1 ^ 2 ~ 1 within a factor of a few. 

The result may be summarized by the simple approximate for- 
mula tt(x) ~ 20(F/0.2) 1/2 a;" 1 , where x = r/R n . This estimate 
is a rather crude approximation (e.g. it neglects the moderate adi- 
abatic acceleration of the collisionally heated jet), yet it demon- 
strates an important feature: tt (x) weakly depends on the param- 
eters of the jet, as long as F S> T„. The estimate of tt also gives a 
simple expression for the photospheric radius, 



R* = t oR n ~ 20 



(e) 



Rn- 



(26) 



The radiation produced by collisional heating in the opaque region 
R n < t < R* is not buried by the optical depth. As demonstrated 
by the radiative transfer simulations in Section 5, it creates a pow- 
erful burst escaping to distant observers. 

4.4 Coulomb heating of thermalized e by ions 

The thermalized e population naturally tends to acquire the so- 
called Compton temperature in the radiation field, Tc, at which 
Compton cooling is balanced by Compton heating due to quantum 
recoil in scattering (e.g. Rybicky & Lightman 1979). If no mech- 
anism heats e , they would quickly reach Compton equilibrium 
with kT e = kTc ~ 1 keV. This however does not happen, because 
thermal e are continually heated by Coulomb collisions with pro- 
tons. As a result, the temperature T e is maintained above Te- 
lls value is calculated below; it satisfies kTc kT e <C m e c 2 in 
the sub-photospheric heating region. 

Nuclear collisions with r rc i > 1 inevitably heat the proton 
component of the jet to a relativistic temperature [j The stirred pro- 
tons acquire a non-Maxwellian distribution with a large fraction 
of protons having kinetic energies <; m p c? . The temperature of 
the accumulated e population is kept at a much smaller value by 
Compton cooling. In this 'two-temperature' plasma, Coulomb col- 
lisions tend to transfer energy from protons to . The thermal ve- 
locity of is well below c and they may be approximated as cold 
(T e = 0) when calculating the Coulomb energy losses of the ener- 
getic protons. A proton with velocity /3 P in the jet rest frame passes 
its energy to the cold e background with rate (e.g. Ginzburg & 
Syrovatskii 1964) 



■ 3. , o- T n±m e c 

£coui = - In A 

2 Pp 



(27) 



where /3 P ~ 1, In A = ln(m e c 2 /hujp) « 20 is the Coulomb loga- 
rithm, and oj p = (47rn±e 2 /m e ) 1//2 . The net rate of energy transfer 
from protons to the thermal e plasma is 



2 

?th ~ ^ In A nCT T n±m e c , 



which gives 



1 



dQth 



nrripC 2 dlnr 



3. .me. 
- ln A — 

2 rrin 



TT W 0.02-TT. 



(28) 



(29) 



It is useful to compare Qth with Q nt h (Section 4.1). From equa- 
tions l|9j, l|29}, and (24) one finds, 



7 Nuclear collisions also create a hot neutron component moving with the 
bulk Lorentz factor V. 
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Y \ 1/2 



(—) 1/2 (—Y 2 (-) 

\5L n J \5T n J \0.lj 



(30) 



Qth 
Qnth 

The thermal and nonthermal heating rates are comparable. 

The e do not keep the heat Qth received from protons. In- 
stead, they immediately pass it to radiation via Compton scattering 
and remain at a temperature kT e <C Qth / n± . The value of T e is 
found from the balance between Coulomb heating and Compton 
cooling of e , 



3 k(T e - T C ) 



2 n± 



Qth, 



(31) 



where tc is given by equation d 1 Ob with 7 ~ 1. We have neglected 
the adiabatic cooling of because its rate is smaller than Coulomb 
heating and Compton cooling rates by the factor tc/icx P <C 1. 
Then we find 



q _ kT e ^ 3m e InA kT c 

e m e c 2 4m p e m e c 2 



0.01 



+ e c , 



(32) 



where O = kTc/m e c? is the dimensionless Compton tempera- 
ture of the radiation field; Be ~ 0.007 for the radiation spectrum 
calculated below (Fig. 5). 



Kompaneets' y-parameter of thermal e is given by 

0.04-tq 
e r 



y = 4T T (6 e - e c 



(33) 



It is comparable to or below unity which shows that Compton cool- 
ing of e occurs in the unsaturated regime. Thermal Comptoniza- 
tion has an important effect on the radiation spectrum, which is 
computed in Section 5. 



4.5 Distribution function of e 
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Figure 4. Momentum distribution of e . Dotted curve shows the case of 
a passively cooling jet (the model with V = 600 from Fig. 1). Solid curve 
shows the case of a collisionally heated jet at the same radius (see text). 
The distribution was calculated at r = 4R n . Vertical dashed line indicates 
the boundary between the thermal and nonthermal parts of the distribution. 
The two parts make comparable contributions to the Compton amplification 
factor A ~ J p 2 (dr^ / dp) dp ~ 1 that measures the average energy boost 
of photons in one scattering by the er^ plasma. The total/integrated optical 
depth tt at this radius is tt = 5; it is strongly dominated by the thermal 
part of the distribution. 



The local distribution function of e is shaped by the processes of 
e injection, cooling and thermalization, and Coulomb heating of 
thermalized e by protons. The distribution is quasi-steady, i.e. it is 
established at a given radius on a timescale much shorter than the 
expansion timescale of the jet. It gradually changes as the jet ex- 
pands. Figure 4 shows the momentum distribution of e at radius 
r = 4i?„ for a typical jet model. The jet has the same parameters 
as in Figure 1 except that it now carries neutrons with F„ = 100 
and L n = 2 X 10 51 erg/s. Synchrotron cooling was neglected in 
this example, i.e. the jet was assumed to be weakly magnetized, 
eb < 1. The distribution has been calculated by the Monte-Carlo 
code described in Appendix B. The temperature of the thermal part 
e ~ 0.03 is consistent with the analytical result (eq. 132b ; it is 
self-regulated so that the balance is maintained between Coulomb 
heating and Compton cooling. The nonthermal part of the distribu- 
tion is formed by the e cascade that results from e injection with 
70 ~ m^/me. 

For comparison, the dotted curve in Figure 4 shows the elec- 
tron distribution that is found at the same radius in the passively 
cooling neutron-free jet. The distribution is Maxwellian, and its 
temperature equals the temperature of the (Planckian) radiation 
field. 



4.6 Radiative efficiency of photospheric emission 

The evolution of radiation density U 1 (measured in the plasma co- 
moving frame) is given by equation 



d\nr V J > \d\nr) 



dQth 
ad d In r 



dQ n 



dlnr 



(34) 



The first term on the right side of this equation describes the adia- 
batic cooling of radiation; it equals — (2/3) C/ 7 in the opaque zone 
and in the transparent zone (the exact behavior of this term near 
photosphere is obtained from the numerical simulation of radiative 
transfer). The second and third terms on the right side represent 
the energy received by e plasma and converted to radiation. Since 
practically all of the energy received by e is passed to radiation, 
these terms effectively serve as sources of radiation energy. Heating 
of the thermalized e population by Coulomb collisions with pro- 
tons dQth/ 'dlnr is given by equation l |29l l. Energy injection into 
the nonthermal e tail dQ nt h/dlnr is given by equation (9j. Sub- 
stituting these expressions to equation l |34t we obtain the equation 
for e = [/ 7 /nrripC 2 , 



de 
dx 



-q(x)e + 



<hh + a nth 



(35) 



where a th = 0.02ro and a nt h = /±r rc i/2 are constants, x = 
r/R n , and q(x) = — (din r7 7 /dlnr) a( j is a dimensionless func- 
tion that equals 2/3 in the optically thick zone tt ^> 1 and ap- 
proaches at the photosphere. The quantity e is the fraction of 
the jet energy carried by radiation; it can also be written in the lab 
frame as 



(36) 



where L 7 = 47rr 2 r 2 t/ 7 c is the isotropic equivalent of radiation 
luminosity, and L = 47rr 2 r 2 nm p c 3 is the isotropic equivalent of 
the jet kinetic luminosity. 
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In the optically thick zone, where q ~ 2/3, equation {35} can 
be solved analytically for e(x), 

R* 



«(*) = 



ei + 3a 

,T 2 /3 



3a 

a; 



1< x < 



Rn ' 



(37) 



where e\ = ej r =jj n and a — a t h + a nt h. The solution may be used 
to estimate e at the photosphere, = R*/R n = To. For jets with 
small ei one obtains 



1 

IV3 



(38) 



r.--= I 0.06r^ /3 + -^i- | | i 

where r ~ 20(y/0.2) 1/2 (Sect. 4.3). Equation {38) estimates the 
net radiative efficiency of collisional heating in jets with r>T„, 
taking into account the adiabatic cooling of radiation until the 
jet expands to transparency^ For example, to = 20 gives w 
0.1 + 0.02r/r n . The two terms represent the contributions from 
the Coulomb heating of thermalized e and the nonthermal e in- 
jection. For typical V /T n ~ 3 — 10 the total radiative efficiency 
= 0.2 — 0.3, with the thermal part comparable to the nonther- 
mal part (cf. also eq.[30}. 



5 RADIATION SPECTRA FROM COLLISIONALLY 
HEATED JETS 

5.1 Thermal and nonthermal Comptonization 

Suppose the jet cools passively at r < R„ and its thermal radiation 
evolves as described in Section 2. The collisional heating begins 
at radius R„ (eq.[6} and quickly loads the jet with energetic e ; 
their typical distribution function is shown in Figure 4. Scattering 
of radiation by the injected e dramatically changes the photon 
spectrum. 

Consider a weakly magnetized jet with Ub "C U-y, when the 
synchrotron cooling of e is negligible. Then the GRB spectrum 
forms via Comptonization of already existing thermal photons ad- 
verted from the center of the explosion. Scattering conserves pho- 
ton number and the average photon energy E in the lab frame can 
be expressed as (cf. eqs.[2]and[3} 



- eYnm p c 
E = 



= eE ^ 4 e r. 



-1/2 r 1/4 



MeV, 



(39) 



where e — L 1 /L is the fraction of jet energy carried by radia- 
tion. The relation d39t is common for all Comptonization models 
of GRBs (e.g. Thompson 1994; Rees & Meszaros 2005; Giannios 
& Spruit 2007). It naturally explains the observed E ~ MeV, as- 
suming a reasonable radiative efficiency e <; 0.1. In the collision- 
ally heated jet, E grows as photons receive energy via two branches 
of heating: thermal and nonthermal (Section 4). The corresponding 
heating rates per photon give 



dE 
d In r 

dE 
dlnr 



0.02td£o — , 
r 



i r p R n 

16 T n r 



(40) 
(41) 



8 Equation 1381 assumes the adiabatic cooling oc r~ 2 / 3 at all r < it*. 
It overestimates the cooling effect - the exact radiative transfer gives less 
cooling (see Section 2). Therefore, equation 1381 underestimates e* by a 
factor ~ 2. 



where Eo is given by equation <[3j . The heating rates and the adi- 
abatic cooling determine the evolution of E(r) in the collision- 
ally heated jet0 However, the known E does not yet determine 
the shape of the radiation spectrum. The spectrum depends on the 
details of Comptonization that need to be calculated. 

The radiation spectra produced by Compton cooled e^ 1 cas- 
cades were previously studied in detail in the context of AGN ac- 
cretion discs (e.g. Svensson 1987; Lightman & Zdziarski 1987; see 
also Appendix A). The model was developed for static sources, and 
Comptonization of radiation in relativistic flows is different for two 
reasons. First, the optical depth evolves as the flow expands. Sec- 
ond, the GRB radiation moves together with the plasma and re- 
mains embedded in it until the jet reaches r ~ Ra ~ F 2 A ~ 



10 lb cm (here A/c 



10 s is the typical duration of GRB 



jets). Collisional heating operates at smaller radii r < Ra, and 
the entire history of heating and Comptonization is 'recorded' in 
the radiation field before it escapes the jet. The spectrum received 
by distant observers is the net result of multi-radius (multi-optical- 
depth) Comptonization in the expanding jet. In this respect, GRBs 
are similar to the relict radiation in the expanding universe. 

The cooling rate of e and their energy distribution at any 
given location depend on the local radiation field. Therefore, the 
evolution of radiation and e plasma must be calculated together. 
This is performed by the numerical code described in Appendix B. 
The code is based on the Monte-Carlo method that solves the ra- 
diative transfer in a jet with self-consistent e distribution function. 
For a given history of heating, the code calculates the evolution of 
e and radiation in the expanding flow and finds the spectrum of 
photons escaping to distant observers. 

As a typical example, consider the jet model from Section 2 
with T = 600, L = 10 52 erg/s, and r = 10 7 cm, but now 
let it contain a neutron component with r n = 100 and L n = 
2 x 10 5 erg/s. The collisional heating in this fiducial model be- 
gins at radius R n ~ 10 11 cm (eq. [6). lust before the onset of 
heating and e creation, the passively cooling jet has E(R n ) w 
1 MeV. The optical depth after the onset of collisional heating is 
T-r(r) = (Rn/r)ro with to ~ 20 (eg. I24t. The heating rates 
in equations d40t and J4U happen to be almost exactly equal: 
(dE/dlnr) nth w {dE/dhir) th w 1.5 {R n /r) MeV. 

Figure 5 shows the spectrum of emitted radiation for the fidu- 
cial model. Although it may not obvious from the figure, the Comp- 
tonized spectrum has two components, which correspond to the two 
parts of the e distribution function (cf. Fig. 4): 

(i) Most photons are multiply scattered by the thermalized 
Coulomb-heated e population and never scattered by the optically 
thin nonthermal tail. This thermal Comptonization dominates the 
emitted spectrum at energies up to 2FfcT e ~ 20 MeV and creates 
the spectrum slope /3 ~ —(2.5 — 3). It corresponds to Kompaneets' 
parameter y ~ 1 that is regulated in the heated jet as discussed in 
Section 4.4. 

(ii) A small fraction of photons are additionally scattered by the 
nonthermal tail of e distribution, which strongly boosts their en- 
ergy. The nonthermal component dominates the radiation spectrum 
at high energies. 

A special feature of collisional heating is that the energy of the 
two spectral components are comparable (eq.|30t. The nonthermal 
component smoothly extends the spectrum through 100 MeV to 

9 Since E/Eq = e for a jet with a conserved photon number, the equation 
for E(r) is immediately obtained from eq. 1351 . 
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Figure 5. Photon spectrum emitted by the collisionally heated jet (solid 
red histogram). The jet has L = 10 52 erg/s, rrj = 10 7 cm, Y = 600 
(same as in Fig. 1), and carries neutrons with T n = 100. Black solid lines 
indicate the slopes that correspond to photon indices a = 0.4 and = 
—2.5. A similar phenomenological spectrum is usually proposed to fit GRB 
observations (Band et al. 2009). The feature near 0.5 GeV is the annihilation 
line. Dotted blue histogram shows the spectrum that would be produced 
if nuclear collisions were 'switched off' at r > 4i? n = 0.2i?,*, i.e. if 
injection was confined to radii R n < r < AR n . The figure does not 
take into account the cosmological redshift of the burst z; the redshifted 
spectrum will peak at(l + z) — 1 MeV instead of 1 MeV. 



the GeV range. It is broad and additionally smoothed by partial 
downscattering in the optically thick plasma before the jet expands 
to transparency. 

Most of collisional heating and Comptonization occurs where 
the jet is still opaque. The Comptonized radiation is released at 
the photosphere i?* and can be called 'photospheric emission' (but 
see Section 5.3 below). The average energy of escaping photons 
in the model shown in Figure 5 is E ~ 2 MeV, which is half of 
Eo ~ 4 MeV. This means that the net radiation efficiency of the 
burst is ~ 50%, i.e. the photospheric emission carries about half of 
the jet energy 

The ratio of the thermal and nonthermal Comptonization com- 
ponents in the observed spectrum is controlled by the parameter 
w w 3t~~ r/r n (the ratio of eqs. 1401 and 141 where to is likely 
to stay around 20 within a factor of a few as long as the condition 
F > T„ is satisfied (Section 4.3). To investigate the sensitivity 
of the model predictions to expected variations in w, we calcu- 
lated three models with equal to — 20 and different T/T n = 3, 
6 and 12. They have w w 0.5, 1 and 2, correspondingly. We found 
similar spectra in all three cases, with slightly different indices 
ft ~ 2.5 ± 0.2 With increasing w, the nonthermal bump becomes 



10 The energy given to photons by collisional heating in the region R. n < 
r < ij* is 2 X 1.5 MeV in the model shown in Fig. 5. Together with the 
initial 1 MeV per photon at R n this would make E = 4 MeV, if there were 
no adiabatic cooling. Adiabatic cooling at r < R+ reduces E by a factor of 
2. 



more pronounced. Large w 2> 1 are not, however, plausible (strong 
nonthermal heating is always accompanied by significant Coulomb 
heating in a realistic jet model). Small w are possible: w can jump 
to zero if r/r n decreases so that nuclear collisions become unable 
to produce pions. This case is discussed in Section 5.4 below. 



5.2 Annihilation line 

The annihilation reaction between thermalized e produces pho- 
tons with energy E' ~ m e c 2 in the rest frame of the jet. The num- 
ber flux of annihilation photons (isotropic equivalent) in the lab 
frame is given by iV ann = 47rr 2 crn ann , where n ann is the density 
of annihilation photons in the jet frame. It obeys the equation, 



47rr ?i a 



rf/Vann 

dr 

Using equations d20t and ( 1241 ) one finds 



dJVann 3tT CTq 2 T 2 



dx 



4 X 2 <7t 



(42) 



(43) 



where x — r/R n > 1 and tq is given by equations J23t . J24t . Inte- 
grating equation | |43I > over x, one finds the net flux of annihilation 
photons emitted to infinity, 

f±Y L 



4r n m e c 2 



(44) 



It is instructive to compare this result with the number flux of orig- 
inal thermal photons in the jet, TV, 



f±Y E 



N 



4T n m e c 2 



2.5x10 



-4 (E^lY 1 ( E ° 

UOOy V0.2J ^MeV 



•(45) 



10~ J cre- 

2 



For our fiducial model shown in Figure 5, Na,nn/N 
ates a rather strong annihilation line that cuts off at E — 2Tm e c 
Most of the annihilation photons are produced well below the pho- 
tosphere. The resulting spectral feature has an extended red wing 
due to Compton downscattering in the sub-photospheric region and 
the variation in the Doppler boost, which depends on the photon 
angle at the emission (or last-scattering) point. 

In strongly magnetized jets, where synchrotron cooling dom- 
inates over Compton cooling, the pair yield Y is reduced (Sec- 
tion 4.3) and the annihilation feature will be weak. 

5.3 7-7 opacity and emission at energies E S> GeV 

To a first approximation, one could neglect the heating at radii 
r S> Rn, and a similar spectrum would be obtained. For in- 
stance, suppose that nuclear collisions occur only in the region 
Rn < r < 4i? n . The result is shown by the dotted curve in Fig- 
ure 5. The spectrum is significantly changed only at high energies: 
the number of photons above the threshold for pair creation is sup- 
pressed. This suppression is caused by the large compactness I at 
small radii, which implies a large optical depth to 7-7 absorption, 

T 77 » 1. 

The extension of the spectrum to ~ 100 GeV in the full model 
(solid curve) is due to the extension of nuclear collisions to large 
radii r ~ 10 3 i?„, where r 77 becomes small and high-energy pho- 
tons are able to escape^The smaller rate of nonthermal heating at 
large r is compensated by the 7-7 transparency at high energies. As 

11 I thank Indrek Vurm for pointing out the effect of continued collisional 
heating at large r on the spectrum shape in the GeV range. 



10 Andrei M. Beloborodov 



a result, ~ 10~ 3 of the jet energy is converted to escaping photons 
with energy comparable to 100 GeV. 

A simple analytical estimate for the optical depth seen by a 
photon of energy E at a radius r is given by, 



,dN/d\nE t 
A-KrcY 2 



(46) 



Here cr 77 ~ 10~ 25 cm 2 is the average cross-section for 7-7 ab- 
sorption by the target photons near the threshold, Et/rn e c 2 ~ 
2r 2 (J5/m e c 2 )~ 1 . This estimate assumes a typical angle 8 ~ 
between the interacting photons. Approximating the spectrum of 
target photons by the Band function with i5 pC ak ~ 1 MeV, one 
finds at E t > -E pca k 



dN 



3 x 10 a 



(A) 

\m e c 2 J 



1+/3 



L 7 ,52 S 



(47) 



din S t 

where L 1 is the total photon luminosity (isotropic equivalent). This 
gives, 



2 x 10 3 



/3-1 



40 

For the typical f3 « 
t 77 = 1) is given by 



9 



7,52 



1 10 GeV J 



1 (JL\ 



2/3 



•(48) 



10 GeV / V600/ 
2.5, the radius of 7-7 transparency (where 



R~n(E) ~ 10 1 



VlO GeV/ 



Goo) 



,52 cm. 



(49) 



.10 GeV, 

The target photons absorbing 10-100-GeV photons have sub-GeV 
energy. Most of them are emitted at r ~ /?* and hence the target 
radiation field at r 2> R* is strongly collimated along the radial 
direction in the jet frame. This creates an "escape cone" for the 
high-energy radiation. The above estimates do not take into account 
this effect. Detailed radiative-transfer calculations (as in Fig. 5) are 
needed to get accurate results. 

A crude estimate for the high-energy luminosity generated at 
large radii may be obtained as follows. The rate of energy injec- 
tion into the nonthermal cascade is given by equation (5}. This 
energy is deposited to radiation via inverse Compton scattering. 
At small radii, the injected nonthermal power is reprocessed by 
the cascade to smaller photon energies E for which r 77 (iJ) < 1. 
As a result, the escaping nonthermal luminosity at energies ~ E, 
dL n th /d In E, is roughly equal to the nonthermal power injected at 
radii ~ R~, y (E), 



1 dL n 



L dlnE 
This gives 

1 dLnth 



J_ _T_ 

16 T\ 



Rn 



a L„ E T 



2-/3 



L dlnE 



(- 



E 



16m„ <r 77 L m e c 2 Y„ \ 2Tm 



1+/3 



(50) 



(51) 



This estimate suggests that at high energies, where 7-7 absorption 
is important, the photon spectrum dN jdE can steepen from the 
slope P to slope /3 — 1. 



5.4 Pure thermal Comptonization by Coulomb-heated 
electrons 

Some GRB jets may have the ratio T/r n near unity. Then nuclear 
collisions are not energetic enough to produce pions and the in- 
jection of relativistic e pairs may not occur. The proton compo- 
nent can still be significantly heated by (elastic) n-p collisions or 
other processes, e.g. internal shocks, and the electron component 
is heated by Coulomb collisions with protons. The electron heating 
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Figure 6. Photon spectrum emitted by a jet with mildly relativistic internal 
motions, which do not lead to pion production. Protons have a mildly rela- 
tivistic temperature in the sub-photospheric region and electrons are heated 
only by Coulomb collisions with protons (see text). 



rate Q t h is proportional to the optical depth tt and can be sig- 
nificant below the photosphere. Radiation in such jets experiences 
pure thermal Comptonization, as the scattering electrons have a 
Maxwellian distribution with temperature T e > Tc (Section 4.4). 

Consider a jet with mildly relativistic protons; they heat the 
electrons according to equation l[29j[3 The resulting radiation 
spectrum is shown in Figure 6. Our simulation assumed the heat- 
ing rate per photon (dE/d In r) t h = 0.08tt MeV and followed 
the evolution of radiation from radius r — 0.05R* where tt = 
20. The average photon energy E = 1 MeV was assumed at 
r = 0.05ii*. It remained close to this value up to r = i?* 
where the jet became transparent and released the Comptonized ra- 
diation. The emitted spectrum is suppressed exponentially above 
E > 2VkT e ~ 20 MeV where kT e ~ 15 keV is the self- 
consistently calculated temperature in the main heating region. The 
example shown in Figure 6 assumes T = 600. We also ran a similar 
simulation for V = 300; it gave a similar spectrum except that the 
cutoff occurred at ~ 10 MeV instead of 20 MeV. 

Any mechanism that keeps protons hot in the sub- 
photospheric region leads to similar Coulomb heating of electrons 
and a similar radiation spectrum. For example, protons may be 
heated by internal shocks in the jet. Internal shocks occur at radii 
r ~ r^in^r, where 8r is the scale of fluctuations and r m i n is 
the Lorentz factor of the slower parts of the jet. The scale Sr 
may be as small as ~ 10 6 cm (the size of the central engine) 
or perhaps even smaller. Then internal shocks begin at a radius 



12 Protons with a smaller, sub-relativistic temperature pass faster their en- 
ergy to electrons (Spitzer 1962). Then, at large tt S> 1, all of the proton 
heat can be taken by the electrons. In this case, the electron heating rate 
simply equals the dissipation rate in the jet. We consider here the case when 
Coulomb collisions are slow enough to create a 'bottleneck' for the heat 
flow from protons to electrons to radiation. In this case, electron heating is 
controlled by the rate of Coulomb collisions. 
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Figure 7. Photon spectrum emitted by a jet whose thermal radiation was 
strongly cooled before collisional re-heating. The jet has the same L, ro 
and T/r„ as in Fig. 5, but its T = 150 instead of 600. 

r ~ 10 10 (r min /100) 2 ((5r/10 6 cm) cm where the jet may have 
a large tt even without production of e ± pairs (cf. eq. |4j V is the 
Lorentz factor of the shocked part of the jet). 

5.5 Strongly cooled and then re-heated radiation 

The standard picture of a passively cooling jet (Section 2) predicts 
that radiation experiences strong adiabatic cooling before reaching 
R n if Rn S> R s , which occurs for modest V. Thus, a regime is 
possible where radiation is strongly cooled before collisional re- 
heating. 

In this case, the emitted spectrum differs from Figure 5, as il- 
lustrated by the simulation shown in Figure 7. It assumes that the jet 
has r = 150, and all other parameters are the same as in Figure 5; 
in particular T/r n = 6, which corresponds to T n = 25. The main 
difference caused by the smaller V is that the radiation temperature 
prior to the onset of heating drops to a low value fcT m i n ~ 60 eV 
(the adiabatically cooled T m i n scales as r 8//3 , see Section 2). It 
corresponds to E m i n « 0.025 MeV and e m i n 6 x 10~ 3 at 
r R n - The collisional re-heating at r <; R n is still strong in 
the model, giving a significant radiative efficiency e ~ 0.4 and a 
significant f«1.5 MeV. The resulting spectrum has a very broad 
peak above ~ 0.1 MeV. It is shaped by thermal Comptonization 
with a large Compton amplification factor A ~ 60. This leads to 
a relatively hard slope f3 ~ —1.4 between -E pca k ~ 0.2 MeV and 
2VkT e ~ 10 MeV where kT e ~ 30 keV is the self-regulated tem- 
perature of the plasma. 

This simulation illustrates an interesting feature of Comp- 
tonization models for GRB emission. Models with efficient re- 
heating do not give the simple Band-type spectrum with the MeV 
break if the mean photon energy E dropped much below MeV prior 
to re-heating. The cooling stage temporarily creates an exponen- 
tial break in the radiation spectrum at a low energy E <C 1 MeV. 
Then re-heating and Comptonization increases E back to ~ MeV; 
however, the recovery of E is achieved mainly by hardening the 



spectrum above the break, with only a minor shift in the break po- 
sition(3 

We conclude that a long stage of adiabatic cooling at r > R s 
has a significant effect on the ultimate spectrum after re-heating. In 
reality, this effect may never occur if the jet does not passively cool 
between R s and R n . Its proton component can be heated e.g. by 
internal shocks. Then Coulomb-heated electrons will keep E and 
e from falling, as discussed in Section 5.4. Thus, it may be that 
even jets with Y ~ 100 — 300 (modest by GRB standards) keep a 
significant e <; 0.1 prior to the onset of inelastic nuclear collisions 
at Rn. Then their emitted spectra will be similar to the spectra of 
high-r jets shown in Figure 5. 

5.6 Impact of synchrotron emission and variability on the 
spectrum 

Two additional effects can change the observed spectra, in particu- 
lar the slope a at E < i5 pca k: 

(i) A strong magnetic field adds synchrotron emission from e 
pairs injected by pion decay. It peaks at energies E s ^ 1 MeV 
(eq. 1 1 8b and can dominate at E < Epc^ because the synchrotron 
spectrum is relatively soft, a — —1/2. 

(ii) Variable jets consist of many thin shells with different param- 
eters, and their radiation spectra vary on timescales as short as 
10~ 4 s (in observer time). The superposition of many different in- 
stantaneous spectra is observed when the true instantaneous spec- 
trum is not time-resolved. This tends to reduce the observed a. 

Models with synchrotron emission and variability will be pub- 
lished elsewhere. 



6 DISCUSSION 

6.1 Formation of GRB spectrum 

Formation of GRB spectrum is a long-standing problem. Much 
of the previous work focused on the optically thin internal-shock 
model (see Bosnjak, Daigne & Dubus 2009 for recent detailed cal- 
culations). The model posits that the observed 7-rays are nonther- 
mal synchrotron emission from electrons accelerated at the shock 
fronts, and introduces phenomenological parameters of this pro- 
cess. Some key issues remained, however, unsettled. Why is the 
nonthermal electron heating efficient? Why do the reported spectra 
of GRBs usually peak near MeV? Why are the low-energy slopes 
of some GRB spectra so hard (a 0, significantly harder than 
e.g. synchrotron emission)? As a possible solution, it was hypoth- 
esized that GRB spectra include a bright photospheric component 
which results from strong sub-photospheric heating (e.g. Rees & 
Meszaros 2005). 

The results of the present paper suggest that the dominant 
component of GRB radiation comes from the photosphere. Col- 
lisional heating naturally gives the photospheric emission a Band- 
type spectrum (e.g. Fig. 5) without invoking unknown parameters 
apart from the Lorentz factors and the initial radius of the jet ro. 
No fine-tuning of these parameters is required to produce the typi- 
cal observed GRB spectrum. The radiative efficiency of collisional 

13 This is a robust result of unsaturated Comptonization in a relativisticalty 
expanding jet. Saturated Comptonization (y 3> 1) would strongly shift the 
peak of the spectrum, however it appears to be not relevant to GRBs as it 
would require photon starvation while realistic GRB jets must advect a large 
number of thermal photons from the central engine. 
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heating is large: it converts i; 30 per cent of the jet energy to es- 
caping radiation. 

In general, collisional heating depends on internal motions in 
the jet. One can imagine three possible regimes: 

(A) The jet is steady (no fluctuations in Lorentz factor) and 
r < r cr it (no neutron decoupling). Suppose also that there is no 
magnetic dissipation. Then the outflowing plasma passively cools 
and emits the quasi-thermal spectrum shown in Figure 1 . 

(B) There are moderate fluctuations ST /T < 1. Then the pro- 
ton component of the jet is heated by internal shocks and nuclear 
collisions between protons and migrating neutrons. Electrons are 
heated by Coulomb collisions with protons, with a well defined rate 
(eq.|29t. The radiative efficiency of such jets is large if the protons 
are hot at radii r ~ 0.1/?* — then Coulomb collisions pass a large 
fraction of proton energy to electrons, and hence to radiation. The 
radiation spectrum emitted by jets with Coulomb-heated electrons 
is shown in Figure 6. 

(C) If there are strong fluctuations 8T/T > 1 or the jet has 
a very high Lorentz factor T > r cr i t , then compound flows form 
with T/r n 3> 1. Nuclear collisions in such jets create pions, which 
leads to injection of e ± pairs with energies ~ m^c 2 ~ 140 MeV. 
The photospheric radius in this regime is regulated by the created 
e pairs. The electron distribution function has an extended non- 
thermal tail, whose shape if determined by the radiative cooling of 
e (Fig. 4). The jet emits the radiation spectrum shown in Figure 5. 
The spectrum extends to very high energies with a slope /3 ~ —2.5 
and has an annihilation line at ~ V MeV whose amplitude is sensi- 
tive to r„. 

Comparison of these theoretical expectations with available 
data suggests that GRB jets are mainly in regime C (and may be 
in regime B in some bursts). The impact of collisional heating on 
the plasma and radiation components of the jet is significant in this 
regime. The jet remains forever dominated by e , n± ~ 20n. The 
produced radiation remains embedded in the jet until it expands 
to r ~ 10 16 cm. Any additional heating processes occur in the 
radiation field already changed by the collisional heating. 

The numerical simulations in this paper were performed for 
weakly magnetized jets, e_g <C 1, where Compton cooling domi- 
nates over synchrotron cooling. Jets with large es are expected to 
have the same photospheric luminosity, with a similar spectrum that 
peaks near 1 MeV but has a smaller low-energy slope a. Numerical 
models for strongly magnetized jets will be published elsewhere. 

We showed numerical examples for a typical GRB jet with 
isotropic equivalent of kinetic luminosity L = 10 52 erg/s. In our 
fiducial model we found R n ~ 10 11 cm, i?* ~ 10 12 cm and 
Rry(E) ~ 10 13 (£/10 GeV) 1 ' 5 cm. The model can be scaled 
to GRBs with different L. Jets with fixed Lorentz factors (e.g. 
T = 600 and T n = 100) and fixed L„/L will have R n oc R> oc 
i? 77 oc L, i.e. the characteristic radii will linearly scale with lu- 
minosity L. The brightest observed GRBs have L ~ 10 54 erg/s, 
which leads to R n ~ 10 13 cm, ~ 10 14 cm, and Rr ll {E) ~ 
10 15 (-E/10 GeV) 1 ' 5 cm. In spite of this big change, the spectrum 
of produced radiation will be similar to that in Figure 5, because 
the ratios Ri,/R n and R^^/Rn are important for the spectrum for- 
mation, not the values of radiiQ The value of Speak is likely to 
increase with L (cf. eq.|39l>. 



The radius of collisional heating is limited by neutron decay at r ~ 
Rp = 3 X 10 15 (r n /100) cm, which does not scale with L. For the bright- 
est jets, the radius of 7-7 transparency at 10 — 100 GeV becomes compa- 
rable to Rp, which could affect the spectrum shape at the high-energy end. 



The slope a of the emitted spectrum is limited by the transfer 
effects discussed in Section 2. The hardest slope found in our ra- 
diative transfer models near 10 keV is about 0.4 (for comparison, 
a Planck spectrum would have a — 1). Practically all observed 
GRBs satisfy this limit (e.g. Preece et al. 2000). However, larger 
q were reported for a few bursts (Crider et al. 1997; Ghiiianda, 
Celotti & Ghisellini 2003; Ryde et al. 2006). This suggests that in 
some bursts the jet may be inhomogeneous on tiny angular scales 

66 < i/r. 

The relativistic jet is causally disconnected on scales 8r > 
?-/r 2 , and different shells 8r can have different radiative history, 
with different Rn and 7?* . Observed GRB light-curves show strong 
variability in a broad range of time-scales beginning from 0. 1 ms. 
The existence of very fast variability is naturally accommodated by 
our model. The shortest timescale of the photospheric emission is 

sum ~ r~ 2 (iVc) ~ io- 4 (r/60o)- 2 i?*,i2 s. 

Note that observations of multi-GeV photons should not be 
used to constrain the radius of prompt emission RmcV as done 
in Abdo et al. (2009a). They derive a minimum V and a mini- 
mum RmcV in GRB 0809 16C assuming that RmcV is the same 
as RgcV (for which 7-7 transparency requires a large value). In 
fact, MeV photons should not be assumed to come from the same 
radius as GeV photons, even when the light-curves in the two bands 
are strongly correlated. The general point is illustrated by the con- 
crete model in the present paper. The same plasma shell that emits 
MeV radiation at R+ can emit GeV photons when the shell expands 
to a larger radius i? 77 (Section 5.3). The temporal correlation be- 
tween MeV and GeV emission is preserved, as photons emitted 
at different radii by the relativistically moving shell arrive to ob- 
server almost simultaneously. There is only a slight delay of the 
very-high-energy component emitted at R 11 . This delay equals the 
observed time of the shell expansion from R* to i? 77 (E), which is 
i? 77 /2r 2 c <> 1 s. 

Similarly, observations of optical radiation that comes from a 
large radius Ro and is correlated with the prompt 7-rays (Racusin 
et al. 2008) cannot be used to constrain the radius of the prompt 
7-ray emission (see also the end of Section 6.3 below). 

6.2 Internal-shock heating 

The jet can be heated in the sub-photospheric region by multiple in- 
ternal shocks as well as by nuclear collisions. A mildly relativistic 
shock front heats protons to a mildly relativistic temperature and 
electrons to an ultra-relativistic temperature if they receive a frac- 
tion e e S> m e /m p of the dissipated energy. At the shock front, 
the electrons acquire a mean Lorentz factor 7i n j ~ e e m p /m e . The 
details of this collisionless process are complicated and can be stud- 
ied by numerical simulations (e.g. Sironi & Spitkovsky 2009). The 
simulations suggest that shocks in a plasma with transverse mag- 
netic field €b > 10 -3 produce a rather narrow (quasi-Maxwellian) 
distribution around 7i n j . 

The volume-averaged e distribution function that results 
from shock heating is similar to that shown in Figure 4. The impul- 
sive heating of electrons to 7 ~ 7i n j is similar to the injection of e 
by nuclear collisions, even though it is concentrated at the propa- 
gating shock front rather than distributed in volume. The heated 
electrons are quickly cooled behind the shock front and create an 
e cascade. 

As a result, the effect of internal-shock heating on radiation is 
in many respects similar to that of collisional heating, and the simu- 
lations in the present paper are useful for understanding this effect. 
If e e > 10 -2 , 7inj exceeds 20. In the sub-photospheric region, 
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where compactness I S> 10, the shape of e distribution func- 
tion at 7 > 20 is not important for Comptonization as the inverse- 
Compton emission from these electrons is anyway absorbed by the 
7-7 reaction. The shape of e ± distribution is quite universal in 
weakly magnetized jets - it is controlled by the development of e 
cascade in the same way for collisional heating and shock heating. 

The thermal part of e ± distribution function must be nearly 
the same in the two cases. It is regulated by the heat supply from 
protons via Coulomb collisions and does not depend on what heats 
the protons - internal shocks or nuclear collisions. Therefore, O e 
is given by equation (32) in either case. 

The main difference between the nuclear collisional heating 
and shock heating is in the rate of electron energy injection. First 
note that the electron energy budget in shock heating is propor- 
tional to e e . A small e e <C 1 implies that the normalization of 
the relativistic tail in the averaged e distribution function is small 
compared to that produced by nuclear collisions, and its contribu- 
tion to Comptonization is smaller. Second, the dependence of shock 
heating on radius is uncertain as it depends on the uncertain vari- 
ability pattern of the central engine. Collisional dissipation has a 
special feature: the 'nonthermal' heating (injection of e by nu- 
clear collisions) dQ nth /dlnr oc r^ 1 scales with r in the same 
way as the 'thermal' (Coulomb) heating dQth/dlnr (Section 4), 
and their constant ratio w is comparable to unity. By contrast, the 
effective w(r) for shock heating may vary, leading to a different 
radiative-transfer solution for the Comptonized spectrum. 

Collisional heating alone gives a 'minimal' emission model. 
Shocks and magnetic dissipation are the usual candidates for ad- 
ditional electron heating, which may create additional components 
of GRB radiation; another mechanism for collisionless heating at 
large radii is outlined in Section 6.3 below. If the nonthermal elec- 
tron population extends to 7 > 100, the relative contribution of 
synchrotron emission increases, as Compton losses are suppressed 
by the Klein-Nishina reduction in scattering cross-section. We did 
not simulate this situation here, and it can be done in the future. The 
general setup of radiative-transfer calculations in this paper and the 
developed numerical code can be used for a broad class of emission 
models - any combination of thermal heating and relativistic elec- 
tron/positron injection in the expanding jet with any magnetization. 

6.3 Heating by neutron decay 

Neutrons carried by GRB jets eventually decay. Their large Lorentz 
factor implies a long decay time F n tp where tp w 900 s. The mean 
radius of /3-decay is 

Rp = ctpT n « 3 x 10 15 (J^j cm. (52) 

The decay occurring inside the jet has a drag effect on the faster 
proton component and reduces its Lorentz factor (Rossi, Be- 
loborodov & Rees 2006). In essence, the decay injects relatively 
slow protons that are picked up by the jet with the relative Lorentz 
factor T rc i = (l/2)(r/r„ + r„/r). This can be described as in- 
elastic sharing of radial momentum between the fast jet and the 
decaying slow neutrons, which decelerates and heats the jet. The 
dissipation efficiency of this process can exceed 100 per cent as the 
jet tends to use the heat to regain its Lorentz factor via adiabatic 
cooling and re-dissipate the energy. 

Most of neutrons decay near the radius Rp. However, a frac- 
tion r/Rp decays at smaller radii r < Rp. Dissipation of ~ 100 
per cent of the jet energy begins at radius Ri ~ (T n /T)Rp. Be- 
tween Ri and Rp, the jet decelerates in the background of decay- 



ing neutrons in a self-similar regime (resembling the deceleration 
of adiabatic blast waves), and its Lorentz factor decreases as r -1 ^ 2 . 
This strong dissipation may generate radiation. 

Note that the decaying neutrons create a perfect maser. The 
new protons injected by /3-decay appear in the plasma frame with 
momentum antiparallel to the flow direction and perpendicular to 
the magnetic field. They immediately begin to gyrate with Lorentz 
factor r ro i and form a ring in the momentum space. This ring is a 
maser that amplifies low-frequency cyclotron waves in the plasma. 
The maser instability develops on a short timescale proportional to 
Wjj 1 where ojb = eB/m e c (e.g. Hoshino & Arons 1991). Damp- 
ing of the excited waves heats the plasma. The waves may also ac- 
celerate particles. Extremely relativistic ion rings were previously 
studied near the termination shocks of pulsar winds and proposed 
to accelerate leptons (Hoshino et al. 1992). A similar heating is ob- 
served in the interaction of comets with the solar wind. In this case, 
a compound flow is formed as the neutral gas around the comet pen- 
etrates the solar wind; ionization of the neutral particles effectively 
injects charges that move with a suprathermal velocity relative to 
the wind plasma and immediately begin Larmor rotation. 

The /3-decay and maser instability produce strong volume 
heating between Ri and Rp. Coincidentally, at comparable radii, 
optical radiation can be released as self-absorption ceases in the op- 
tical band. Besides, the jet becomes transparent to very high-energy 
photons. Thus, interesting radiative signatures may be expected. 
The radiative efficiency is, however, uncertain and likely smaller 
than the photospheric e* ~ 0.3 — 0.5. The emission will occur in 
the optically thin regime and can be of the type modeled by Stern 
& Poutanen (2004) and recently by Vurm & Poutanen (2009). The 
study of possible radiative signatures of neutron decay between Ri 
and Rp is deferred to a future work. 

One feature of emission powered by neutron decay can be 
predicted. The emission will arrive to distant observers with a 
slight delay with respect to the photospheric emission produced 
by the same neutrons via the collisional mechanism at r ^ i?*. 
As the jet expands from i?* to r S> R*, a neutron with Lorentz 
factor T n shifts with respect to the plasma jet a radial distance 
Ar « (r - R*)(r- 2 - r~ 2 )/2 « r/2r 2 , which corresponds 
to observed delay 

This neutron-drift delay appears to be consistent with the detected 
delay A£ D b s ~ 2 s of the prompt optical emission with respect to 
the main GRB pulses in the 'naked-eye' GRB 0803 19B (Beskin et 
al. 2009). 

The delay of the GeV source detected by Fermi can have a 
similar, geometrical reason: if it operates at radii much larger than 
the source of MeV emission, its emission must be delayed. 
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APPENDIX A: ELECTRON-POSITRON CASCADE 

The compactness parameter I is large in the sub-photospheric region r < 
Ri, where most of collisional heating takes place (Section 4). It implies 
fast Compton cooling of er^ and quick 7-7 absorption of energetic pho- 
tons. Therefore, the timescale for the cascade development following the 
injection of an electron (or positron) with 70 ~ m w / m e is short compared 
with the jet expansion timescale. Then the distribution function of nonther- 
mal e^, dn±/d*y, is locally (at a given r) qausi-steady and satisfies the 
equation 



d_ / dn±\ 
d^y \ d-y J 

where S(-y) = 



5(7), 



(Al) 



dh±/d r y is the creation rate of secondary e , and 
m e c 2 7(7) is the energy loss rate of electron (or positron) with a Lorentz 



factor 7. The energy loss is due to Compton scattering, synchrotron emis- 
sion, and Coulomb collisions with thermal electrons and positronsF^I 
Coulomb collisions dominate at small 7 fa 1 (and lead to quick fhermaliza- 
tion of non-relativistic e*). In this Appendix, we focus on the relativistic 
tail of the distribution, where Coulomb collisions are negligible com- 
pared with Compton and synchrotron losses. Then, 



m e c 2 7 a -_ (JTC ([/ KN + U B )(T - 1), 



(A2) 



where Ckn(7) is the energy density of photons with energy E' is 
m e c 2 /7, i.e. below the Klein-Nishina cutoff in the scattering cross section. 
The source function 5(7) may be written as 



S(7) = 9irf)Moh, M 



f±r Icl m p 

70 me 



(A3) 



Here h is the rate of nuclear collisions in the compound flow (eq. [7), and 
Mo is the multiplicity of primary injected with Lorentz factor 70 fol- 
lowing a nuclear collision. Then 9(7) is a dimensionless function that rep- 
resents the source of secondary created by one primary e~ or e + . This 
function is calculated numerically using Monte-Carlo simulations of the 
cascade. 

Integration of equation jAlt yields 



^T nth _ <7 T r dn± 

r 



d-y 



dr/ 



3 kJ nniec 2 G( 7 ) 

8 M0Tn U KN + U B —l' 



G(7) 



9(7') drf- 



(A4) 



(A5) 



When Ub "C U-y, this equation becomes 



Gil) 



(T v r dn+ 3 n ... 

~~T^ = o -^0 T n v 2 " , (A6) 

r d-y 8 n 7 e(7)(7 — 1) 

where n 7 /n ~ 10 5 is the photon-to-baryon ratio (the main parameter of 
the GRB jet, see eq.[2), and 5(7) = UKN/ n j m eC 2 represents the mean 
dimensionless energy per photon below the Klein-Nishina cutoff. The typ- 
ical e in the calculated models is near 3 X 10~ 3 and varies slowly with 
7. The dimensionless function G(-y) equals the number of secondary 
injected with Lorentz factor above a given 7 in the cascade triggered by 
one primary particle. In particular, G(l) = Ms is the total number of sec- 
ondary e . The function G( 7 ) decreases from G(l) = Ms to 0(7) ~ 1 
at 7 iS 70 , which implies a relatively slow dependence on 7, with the aver- 
age slope din G/d In 7 — \nMs/ In 70 ~ —0.7. Then equation iA6\ 
implies that the optical depth of the nonthermal population r n th sharply 
peaks at 7 ~ 1. Its value is small, T nt h <C 1 (and much smaller than the 
optical depth of the thermalized e^). However, the effect of nonthermal 
population on radiation is measured not by r nt ^ , but by the Compton am- 
plification factor A = J 7 2 dT nt i 1 . The amplification factor peaks at large 
7. 



dA 3 n 7 G( 7 ) 

— = -M t„ — -, 

a In 7 8 n 7 £(7) 



(A7) 



Particles with 7 e 2> 1 generate photons that are absorbed by the 7 — 7 
reaction. As a result, in the main heating region r < R*, particles with 
7 ii 20 contribute to the development of the e*- cascade, while particles 
with 7 Si 20 shape the scattered radiation spectrum. 



APPENDIX B: NUMERICAL CODE 

The code is designed to simulate the self-consistent evolution of the radi- 
ation field and the e plasma in the jet. For collisionally heated jets con- 
sidered in this paper, the following quantities are known at all radii: the 
injection rate of primary with 70 ~ 300 (eq. |7J and the corresponding 
energy injection rate (eq.|9j, the density of the accumulated thermalized eP 1 



5 We neglect in this paper the possibility of energy exchange between 
thermal plasma and nonthermal due to collective processes. 
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component (eq. I24t and the heating rate of this component (eq. 1291 . The 
code aims to find the temperature of the thermalized v. : population T e (r), 
the nonthermal tail of distribution, and the radiation field at all radii. 

The calculation is split into two parts: (i) global radiative transfer in 
a jet with a given distribution function, and (ii) calculation of dis- 
tribution function for a given radiation field. The consistency between parts 
(i) and (ii) is reached via iterations as explained below. Note that part (i) is 
a global problem, while part (ii) is local and can be solved separately at all 
radii. Temperature T e and the nonthermal tail of er^ distribution at a given 
r are determined by the local radiation field, Coulomb heating rate, and 
injection rate. 

Radiation has a Planck spectrum at early stages of jet expansion (i.e. 
at small radii), with the temperature determined by the initial size of the jet 
and its energy. In the simulations, the initial thermal radiation is sampled by 
a large number of Planck photons (~ 10 9 ), which are injected at a small ra- 
dius and their scattering is followed until the photons escape. The code can 
also simulate the injection of synchrotron photons and follow them together 
with the Planck photons. In this paper, only weakly magnetized simulations 
are presented, where Compton scattering is the dominant mechanism of 
spectrum formation, and synchrotron emission is neglected. Before the jet 
expands to transparency, the photons are multiply scattered and may be ab- 
sorbed by the 7-7 reaction. In each scattering event, the scattering electron 
is randomly drawn from the local distribution function, and the exact 
Compton cross section is used to randomly perform the scattering. 

The radiative transfer is calculated in the static lab frame, assuming 
that the plasma flows in the radial direction with a bulk Lorentz factor T. 
Since T is large (10 2 — 10 3 in the simulations) essentially all photons flow 
outward, and most of them have tiny angles 9 ~ T -1 with respect to the ra- 
dial direction. The radiation is essentially comoving with the plasma flow. 
Therefore, one can view the transfer as the evolution of radiation in time 
t = r/c — time and radius are almost equivalent choices for the inde- 
pendent variable in the problem. Between successive scatterings at radii ri 
and T2, the photon propagates along a straight line in the lab frame, and 
its angle with respect to the local direction of the radial jet, 9, changes: 
sin 62 = ixi/rz) sin Si. This change automatically (and exactly) de- 
scribes the adiabatic cooling of radiation in the opaque zoneF^I 

When solving the radiative transfer with a trial T e (r), we find the en- 
ergy gained by radiation (per photon) from scattering on thermal e ■ This 
is done by defining a radial grid r ; and accumulating statistics of scattering 
in each bin A In r during the Monte-Carlo simulation of the radiative trans- 
fer. Thus, we evaluate (dE/d In r)th( r ) for our trial model. If it exceeds 
the required (dE/d In r) t h (given in eq. 1401 we reduce T e (r) in the next 
iteration. 

The nonthermal tail is given by equation (AlJ, which contains the 
source function 5(7) with shape 3(7) (eq, IA3t . We find 9(7) numerically 
using the Monte-Carlo simulation of the cascade in the local radiation field 
(which is known after calculating the radiative transfer in the previous itera- 
tion). The distribution of nonthermal at small 7 is affected by Coulomb 
collisions with the thermalized population. This effect is included by 
adding the Coulomb losses to 7 in equation )AU . The losses are evaluated 
approximately by assuming a cold background!* 7 ! they are given by 



equation 1271 (Coulomb losses are similar for energetic protons and e^). 
The updated nonthermal tail and T e (r) are used in the calculation of ra- 
diative transfer in the next iteration. 5-10 iterations are usually sufficient to 
accurately find the self-consistent solution for the radiative transfer and the 
e distribution function. 

The iterative method also allows the code to achieve a self-consistent 
treatment of 7-7 absorption, which is a nonlinear effect. When the code 
solves radiative transfer, the opacity to 7-7 absorption k 77 is evaluated us- 
ing the radiation field saved from the previous iteration (or an initial guess, 
for the first trial). In each radial bin A In r, the radiation field is saved on a 
grid in the (9, i?)-space as a collection of ng xnj; 'monocromatic beams' . 
The 7-7 opacity seen by a given photon with energy Eo propagating at 
angle 9q is calculated by integrating over all target 'beams' (A9, AE) 
(with a random azimuthal angle) that are above the threshold for reaction 
7 + 7— > e+ + e _ . As the photon (Eq,9o) traverses the radial bin A In r, 
the probability of its survival is cxp(— k 77 Ar / cos 0q). 



16 This can be understood by considering the toy problem of coherent and 
isotropic scattering in a cold jet. Then a scattering event does not change 
the photon energy in the local plasma frame E'; it only changes its angle. 
Between successive scatterings, the energy of the freely propagating photon 
in the lab frame E = const, and E' = EF(1 — /3 cos 9) is decreasing be- 
cause of decreasing 9. In addition, the propagating photon becomes prefer- 
entially beamed outward in the plasma frame (9 1 decreases). Next scattering 
again randomizes cos 9' and destroys the preferential beaming, suddenly in- 
creasing (on average) 9' . As a result, the next scattering on average reduces 
the photon energy E = ET(1 + /3 cos 9') in the lab frame. 

17 The exact shape of the distribution function in the region connecting 
the thermal and nonthermal parts requires the full treatment of Coulomb 
collisions with a finite-temperature plasma. However, this region radiates 
very little, and the approximate matching of the thermal and nonthermal 
components (as in Fig. 4) is sufficient for the radiative transfer simulations. 



